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HYDROMAGNETIC WAVES AND DISCONTINUITIES 


IN THE SOLAR WIND 


I. Introduction 

Studies of the interplanetary medium have produced results of astrophysical and 
geophysical significance. During the past ten years it has also become possible 
to treat the interplanetary medium as a laboratory sample of a cosmic plasma. 
One of the major results of these experimental investigations is the demonstra- 
tion that hydromagnetic theory is applicable to this tenuous plasma. Of course, 
fluid processes represent only a subset of the great wealth of phenomena occurr- 
ing -and one must remember that the applicability of hydromagnetic theory is 
limited. Moreover, these limitations are presently the subject of great interest, 
and the transport parameters remain to be determined by the development of an 
appropriate kinetic theory. But the applicability of hydromagnetic theory to the 
solar wind is very extensive nevertheless. While the current trend in the liter- 
ature emphasizes the limitations of the fluid theory, this review emphasizes its 
positive results. 

Various modifications of the hydromagnetic theory for a single species fluid 
having an isotropic pressure have been made for specific solar wind problems. 
We shall develop a general form of hydromagnetic theory appropriate for the 
solar wind, and we shall show how the various special theoretical results can 
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be unified and systematically derived from this general formulation. We shall 

* 

also show how solar wind plasma and magnetic field observations support and 
can be* organized by this theory.' The emphasis Is thus on the basic theoretical 
and experimental results*. and relations between them. Tor this- reason,, complex, 
non-reproducible phenomena and descriptive results very specific to* the solar 
wind are not discussed. Likewise, theoretical work which has no present con- 
nection with observations is not reviewed. 

II. Summary of Solar Wind Properties. 

Since there exist g*o,od recent reviews of the- observation material-,- (Axford; • 

1968; Hundhausen, 1968, 1970; Ogilvie, 1970; Scarf, 1970; Vasyliunas,. 1971) we 
shall not go into great detail in the description of specific experiments and their 
results-. The aim of this- section is to provide the reader- with -a- condensed' sum- . 
mary of those properties presently established together with .the necessary , 
references to the papers which give the detailed information not- presented* 
here. 

An exhaustive review of experimental methods, drawing attention to the im- 
portant conversion of engineering quantities (count rates, currents, etc.) to 
physical quantities (temperature, densities, etc.) characteristic of the particle 
distribution functions has been published by Vasyliunas (1971). As- he points out, 

k b am mu si t miete ti i§ ite d&timmtiea eltfe Mgffe sf $8 
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present in various parts of velocity and configuration space. The results, 
though always an approximation, represent the essential features of f(r , v , t), 
defined such that f d r dv is the probability of finding a particle in d r dv at 
time t. 

In the fluid description of the plasma, it is considered as a conducting fluid, 
whose motion and properties (for example transport properties) are modified 
by the presence of electric and magnetic fields. The state of this. fluid at a 
given time is specified by the parameters, density, bulk speed etc. The evalu- 
ation of these quantities in terms of the distribution function is covered in many 
textbooks (for example Boyd and Sanderson, 1969), where it is shown that 


n (r , t) 


f (r , V, t) d V 


( 2 . 1 ) 


u (r , 



f d v 


( 2 . 2 ) 


TO. 


p 

3 n k 


(2.3) 


where P is the pressure tensor, defined in terms of the second velocity moment 
of f. 

We now discuss the observed properties of the plasma particles in terms of 
these quantities. For this purpose it is important to introduce the idea of scales 
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of length, (or time). The idea of "macroscale," "mesoscale" and "microscale," 
introduced by Burlaga and Ness (1968) and discussed further by Burlaga (1969) 
has been successfully used to classify and interpret data, a length being con- 
vertible to the corresponding time by dividing by a representative bulk speed. 
The mesoscale deals with lengths of £1AU and the microscale with lengths of 
£0.01 AU respectively (see Figure 1). It is clear that mesoscale changes in the 
fluid parameters are indicative of changes in the flow field on a large scale. 
Microscale changes reflect the passage past a given observatory of waves, 
discontinuities or rapid changes. Examples of such effects will be discussed 
later at length. 

Observations of fluid quantities on the macroscale show that while there is 
considerable variability, the bulk speed is seldom below 300 km sec"i and the 
low speeds are between 300 and 350 km sec -1 . At these tim.es the temperature 
is about 3 x 10 4 °K, so the flow is supersonic, with Mach, num.ber of 
order ten. For these conditions, referred to as the "quiet solar wind" and 
usually used to compare with theoretical predictions at 1AU, the density is = 8 
cm -3 . At other heliocentric distances, between say 0.25 AU and 5-10 AU, the 
flow speed is predicted to remain closely constant, the density decreases as 
l/r 2 ,and the temperature decreases as r increases at a rate whose details de- 
pend upon the heat conductivity of the plasma. It is expected to be somewhat 
slower than the adiabatic rate. Relatively little observational material 
has been presented about conditions at heliocentric distances other than 1 AU. 
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The proton temperature is anisotropic about the direction of the interplanetary 
magnetic field, T )( /T x being 2:1, which is less than would be expected -from the 
Coulomb collision length in the medium alone. 

The electron density is such as to preserve charge neutrality, and electrons 

convected at the same bulk speed as the protons. The measured electron 

temperature is essentially constant at a value of 1 to 1.5 x 10 S °K, with an ani- 

•*1 

sotropy about the field direction of ,<1.2. This means that the electron distribu- 
tion function is quasi-stationary in time, and that the flow of the electron gas is 
subsonic, since the electron thermal speed is several thousand km/sec. 

At the present time the density, temperature and its anisotropy, the bulk speed 
and heat flux have all been established for both protons and electrons at 1AU. 

These values for the quiet solar wind as defined above are given in Table -I. The 
values of some derived quantities are also given. 

The value of /3 s nk(T + T )/(B 2 /Sn), in the solar wind at IAU is typically close 
to unity. Other species of particles exist in the solar wind. Extensive observations 
show a highly variable admixture of helium. If one disregards the short periods when 

the proportion of helium to hydrogen exceeds 10%, which seem to be associated 
with solar activity, we see that one can generally neglect the presence of this 
second species . The ratio of the helium thermal energy density to that of the 
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Table I 


Observed and Calculated Properties in the Quiet Solar Wind 


Observed 


Plow speed 

320 km sec -1 

Proton and Electron density 

8 cm -3 

Proton Temperature 

4 x 10 4 °K 

Proton Thermal Anisotropy Ratio 

1.9 

Electron Temperature 

1-1.5 x 10 s °K 

Electron Thermal Anisotropy Ratio 

1.1 

Magnetic Field 

5y 

Proportion of Helium (n a /n p ) 

0.05 

Heat Conduction Flux (total) 

0.01 erg cm' 2 sec" 1 

Calculated 


Proton flux 

2 x 10 8 cm" 2 sec" 1 

Kinetic Energy Flux 

0.22 erg cm" 2 sec -1 

Kinetic Energy density 

7 x 10" 10 erg cm -3 

Proton Thermal Energy Density 

6 x 10" 11 erg cm -3 

Electron Thermal Energy Density 

1.5 x 10" 10 erg cm -3 

Magnetic Field Energy Density 

10 " 10 erg cm" 3 

B 

1 

— 0.3 


-1.0 

Variability at 1 AU 


Flow Speed 

275-800 km/sec _1 

Density 

~ 1.0 to ~ 50 cm" 3 

Proton Temperature 

< 3 x 10 4 to 5 x 10 5 °K 

Proton Thermal Anisotropy Ratio 

1 to 3.0 

Electron Temperature 

1 to 2 x 10 5 °K 

Proportion of He (n a /n p ) 

0.01 to 0.20 

Helium/HydrogenBulk Speeds 

<5% 
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hydrogen is 0.2; that of the two kinetic energy fluxes is ;0;2,->but since the bulk 

t 

speeds seem always to be equal within the accuracy of observation (a few per- 
cent) the principal effect of the helium is to increase the fluid .particle density 
and pressure. The possible effects of the admixture of a second species on the 
stability of the plasma seem not to be important in the fluid approximation, 
{Fredericks and Scarf, 1965), -The other species recently detected in the solar 
wind at quiet times are present in such small amounts that their effect is prob- 
ably negligible. 

(Bame et. al., 1970). We are dealing with a multi- component, anisotropic, colli- 
sionless plasma, with supersonic protons and subsonic electrons, which has 1. 

HI. Basic Equations 

A. Discussion. As shown above, the general equations 'for fluid motions in 

the solar wind must include the magnetic field, the particle anisotropies; and the 

different species. There are 5 transport equations (conservation of mass (1), 

momentum (3), and energy (1)) for each species and those are coupled with one 

another and with the equations for the current and magnetic field. Because of 

the complexity of this system, it has not been studied in its general form. Two 

approximations are commonly used: -(1) a single fluid approximation, obtained 

by averaging over particle species, which is ^ valid for an anisotropic, hydromag- 

netic fluid; and a 2-fluid model for an isotropic fluid with.B = 0. , The 2-fluid. 

model was developed for investigations of macroscale solar wind properties. 

< ' 

The results of this paper are related to the single fluid approximation. 
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B. Single- Fluid Approximation 


1. Conservation of mass. Averaging the mass conservation equations 
for each species gives 


+ V* ( p v) = 0, 

0 t 


where 


and 



(3.1) 


(3.2) 


(3.3) 


here is the particle mass, n 1 is the density of particles of species i given by 
(2.1), and u. is the average speed of species i given by (2.2). 


2. Conservation of momentum. For the solar wind, the momentum 
equation is 

~ (p y) + V- (P + P em + p v v) = 0 (3*4) 


where P em is Maxwell's stress tensor. 


P 

~ em 


5... JL-hh 

1 1 8 77 4 77 


(3.5) 


and f is the mechanical stress tensor summed over particle species, which can 
be written in the form 




B. B. 


^kll “ PH ) + Pt-f 3, ; 


kL 


B 2 


kl ij 


(3.6) 
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where p is the pressure of the kth species perpendicular to B, 


p ki = n k m k 


w ki w kL f d v 


(3.7) 


and p, .. is the pressure parallel to B, 

“ 13 


p kii 


n,. m. 



w kll 


f dv; 


(3.8) 


where 


w k = 


- u. 


(3.9) 


One might add a viscosity term to P, but it seems to be negligible for macro- 
scale processes in the solar wind. Putting (3.5) and (3.6) into (3.4) gives 


x V . 
3t 


O v V 4- I 



= 0 


(3.10) 


where 




YL (Pk " ~ 

k 

(B 2 / 4 77) 


Noting that 


B (p v) 

3 t 


+ V' (p V v) = 



and using the condition 7 • B =0, one gets from (3.10) 


(3.11) 
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(3.12) 


(B • V)B ^ 

4 77 

if 4 = constant. This equation with £ = 1 was first derived by Parker (1957) 
from particle orbit theory. Clearly, g is a measure of the anisotropy, being 1 for 
an isotropic plasma and <1 for an anisotropic plasma with p > p . Note that 
the anisotropy enters only in the (B • V )b term. Since (jj • V )b is proportional 
to the curvature of the field lines, it is clear that the anisotropy directly affects 
the motions if and only if the field lines bend. 

Equations (3.10) and (3.12) are fluid equations, and they imply a collective inter- 
action of the fluid particles. For motions normal to B this is caused by the 
magnetic field; collisions are not necessary. This was shown by Burgers (1960), 
Chandrasekhar (1960), and Whang (1970) , who obtained (3. 10) for for a collision- 

less plasma. For motions along B there must be collective plasma interactions 
if the motions are to be fluid-like. These motions are not yet hilly understood. 

It is possible that on some scales one should use a kinetic equation for motions 

r 

along, B , rather than (3.10). 

• j 

3. Frozen- Field Condition. Assuming infinite electrical conductivity 
and neglecting the effects at ion and electron inertia, one has E = - v x B/c , 
where c is the speed of light.. In this case v and B are related by the equation 

.uV x (-y * B} (3.13) 
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This is essentially the same as a form of the Beltrami equation for vorticity in 
a non-viseous, incompressible medium (Serrin, 195 9, page 151), and it implies-- 
that the lines of force (or vortex lines in the case of the Beltrami equation) are 
"frozen" to the fluid. It was first discussed by Alfven (1943). 

4. Closure. For many purposes, it is not necessary to consider higher 
moments of the distribution function. However, the conservation of mass and 
momentum equations ((3.1) and (3.12)) and the frozen-field equation, (3.13), give 
only 7 equations whereas there are 9 unknowns (v , B , p , and p |( ) . Thus, 2 
additional equations are needed to obtain a closed set of equations describing the 
plasma motions. 

For macroscale studies, it is customary to assume isotropy ( £ = 1), which 
reduces the number of variables to 8. Then only one additional equation is 
needed, and it is customary to use the adiabatic approximation p = A(S) y which 
gives a relation between the pressure and density which is valid for constant 
entropy, S. In the kinetic theory of gases it is shown that y - (N + 2)/N where 
N is the number of degrees of freedom. The application of this equation to the 
solar wind has no theoretical or experimental justification, but it is obviously 
mathematically convenient. 

Another means of closing the equations, which is sometime applied to the solar 
wind when discussing problems involving anisotropy, is the method of Chew, 
Goldberger and Low (1956). Their closure equations are 
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(3.14) 


and 




d /Pii b2 \ 

dt U ) 


= 0 


(3.15) 


These equations assume no heat flux, but can be extended for a constant heat 
flux. The CGL closure equations were derived on the assumption that ft = n k T/ 
<B 2 /8 t t) « 1. This assumption is not applicable to the solar wind. 


We shall close the system with the following equations 

Phi = V pj* (3.16) 

Phil = A. ii pI‘ < 3 - 17 > 

The appropriate values for y L and are unknown. An obvious choice for 
microscale processes where the field is relatively orderly is y y = 2 (because 
there are 2 degrees of freedom normal to B ) and = 3 (because there is 
obviously just one degree of freedom along B ). We cannot exclude the possi- 
bility that compressions along B are isothermal since the conductivity along B is 
high, especially for electrons .(p « p)\ in this case = 1 rather than 3. One 
should also consider the possibility that y depends on the scale, since on a large 
scale the field is more disordered (perhaps turbulent) than on a small scale. 
These are problems for future work. 
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5. Energy Equation. The forms of energy to be considered to drive the 


microscale processes in the solar wind are the kinetic or streaming energy, the 
parallel and perpendicular thermal energy and the magnetic energy. Let 


The energy flux is 





(3.18) 


F = 



+ P 


v + S -f q 


where S is the Poynting flux 


(3.19) 


S 


(v x B) x 8 
4 77 


c 


(3.20) 


and q is the heat flux 

d v w? Wj f. (3.21) 

The heat flux, q , has been neglected in most of the dynamical studies of micro- 
scale processes, but this assumption has not been justified experimentally and 
there are recent studies which indicate that electron heat fluxes may be very 
significant under some circumstances (Burlaga et al., 1970, Hundhausen and 
Montgomery, 1970). 



13 



The conservation of energy is expressed by the equation 


5-l + V-F = ° (3.22) 

6. Remarks on Velocities. It is important to distinguish the various 
velocities which are used above. The basic equations involve the center of mass 
velocity, v , which is a momentum average over the particle -species, viz., as 
pv = 2 a tn.u i . Here p is the average density, p = 2 re m i5 and u i is the average 
flow velocity of the ith species, 

n C 

U. = _L V. f dv- (3.23) 

1 n 1 1 

V 


where v . is a point in the velocity space of the ith species. The momentum 
tensor, given by (3.6), and the heat conduction vector, given by (-3 .21), involve 
the thermal velocity measured relative to the center of mass velocity, i.e., 
w. = v. - v . This presents a problem since in practice one would usually 
determine P and q for each species separated, by integrating over the thermal 
speed w' = , measured relative to v j . Since particles of all species 

are frozen to the lines of force, they should all move normal to B with the same 
velocity, i.e., v = v , but there seems to be no reason why different species 

■*- i 1 

could not move parallel to B with different speeds. However, observations 
(Ogilive et al., (1968); Robbins et al., 1970) indicates that in the solar wind the «’s 
and protons have essentially the same mean speed (within 10 km/sec), and the 
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nean speed of the electrons does not significantly effect v -because of the small 
electron to proton, mass ratio, so for most purposes we can set v = u i for solar 
vind particles, in which case w.' =« w i . 

7. WAVES 

A. Introduction. Although the theory for hydromagnetic waves in an iso- 
;ropic fluid is given in many texts, the theory for an anisotropic, multifluid, 
jollisionless plasma is not well known, and differs from the isotropic theory in 
some important respects. Since the anisotropic theory is the relevant one for the 
;olar wind, we develop it in some detail. We shall consider a number of special 
;ases which show the basic physical properties of the waves. The general 
situation is discussed briefly at the end of this section. 

Ine can distinguish 3 classes of waves: 


1. 

v B i 0, 

v l = 0 

(4.1) 

2. 

v li = °’ 

v^O 

(4.2) 

3- 

• v ll * 0 ’ 

• V|^ d 0 • 

- ' (4.3) 


Here, and in what follows in this section, v is the perturbation velocity due to 
;he wave. ' - ' * - ' ' ’ 

B. Waves with vjj £ 0, vj. =0. hi this case, the oscillations are only along 

" “ 7 7 : ; * . . ; /*» 

;he magnetic field. Equation (3.13) shows that B B/Bt = 0 for such waves, i.e., 
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they involve no variations B , Since B = constant and is also constant, the 
linearized form of (3.10) gives 



(4.4) 


where c jj s 3p k|| /B is the speed of sound of the kth species. If p k = A p k y, 
then c 2 = y k T k /m k . Let 

* f 

c i ^p~ x YL (c ^i (4 * 5) 

k 

For a plasma consisting of electrons and protons, 



y + TJ +,y_ TT 


nr 


(4.6) 


If Cm is a constant, then (4.4) and the continuity equation give 


B 2 


v„ 


- c fj V„ 2 v„ = 0 


(4.7) 


This is an equation for sound waves propagating along B with speed c 2 . The 

r\j II 

wave is illustrated in Figure 2a. 


The existence of acoustic waves in a eollisiohless plasma is somewhat para- 
doxical. However, such waves (ion acoustic waves) are predicted by collision- 

• » * 

less theory and are expected to propagate with little attenuation when T e > > T. . 
The restoring force is provided by am electrostatic effect. 
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Such waves have not been identified in the solar wind, and their coupling with 


other modes (assuming ion acoustic waves do propagate) is poorly understood. 
Such waves do not appear in Whang's model for flow past the moon since he 
assumes that particle motions along B are governed by orbit theory rather than 
by the flow equations. 

C. Waves with Vj_ = 0, v M / 0. The equation of motion (3.12) can be written 


where 


and 



U = - 


4 77 


(B -V) B. 


(4.8) 

(4.9) 
(4.10) 


Clearly, G is due to gradients in the total perpendicular pressure. M is caused 
by the tension in the magnetic field lines, which is non-zero if and only if the 
field lines are curved. Note that the anisotropy appears only in M. We can 
distinguish 3 cases, shown in Table H. (Bear in mind that we are assuming v M = 0. 

Iso note that we are using the terms magnetoacoustic and fast wave in a limited 
sense.) Now, consider each of these cases in turn. 
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Table II 



G 

M 

Alfven wave 

0 

7^ 0 

Magneto acoustic wave 

^ 0 

0 

Magneto "fast wave" 

^0 

4 o 
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1. Alfren wave (G = 0, M £ 0; Vj| = .0). In this ease .there .are no pres- 
sure gradients so n, T, and B are constant, but the magnetic field lines oscillate 

♦ /-o ’ * ^ < 

like waves on a string. Let us choose a coordinate system such that 9 is along 

~ ■ •* 

^ and take the perturbations along y : S B a b (z) y , S v = v(z) y (see Figure 3) . 

Consider the case when B is parallel to k. The force equation (3.12) gives for a 

homogeneous plasma 


_ I - _8_\ B_b 

3t \ 4 7 t/3z 


(4.11) 


This implies that §B is parallel to 8 v , The frozen- field condition, gives 


3b_g3v 
d t d z 


(4.12) 


Combining these two equations gives the wave equation 


3 2 b _ / B 2 A B 2 b ^ 3 2 b 

a! 2 \4 77 P j = V A 


(4.13) 


Thus, the waves move with the speed along b . Assuming perturbations of the 
form b 0 exp [i(kz - cj t)] , we obtain the dispersion equation 


&,2 _ v' 2 k 2 = 0 (4.14) 

A 

Note that the speed of the waves is affected by the anisotropy; the larger the 
anisotropy, the slower the wave, all other things being equal. Clearly, the waves 
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' will propagate only if V A is positive. For an isotropic plasma, g = 1 and V A = 
’ V A , which is the familiar result from magneto hydrodynamic theory. 


From the frozen-field equation (3.13) 

U> _ _ v , AA 

B t a ^ z 


(4.15) 


which with (4.12) gives 



This can be written. 


(4.16) 


Unti and Neugebauer (196 ) searched for sinusoidal Alfven waves in the solar 
wind using magnetic field and plasma data from Mariner II. One result is given 
in Figure 4. This shows B (smoothed by a 19 point running average) in a co- 
ordinate system for which the minimum field fluctuation was in the xy plane and 
the average field along z, the bulk speed of the plasma is also shown. One sees 
a nearly sinusoidal, plane polarized wave with B B strongly correlated with V 
and § B normal to B . Using the measured B and n to compute 8 v they found good 
agreement with the measured bulk speed perturbations (Figure 5), suggesting 
that (4.17) is satisfied. The initial density measurements, computed on the 
assumption of a strictly radial bulk velocity, did not satisfy the condition n = 
constant, necessary for an Alfven error, (in fact the density was strongly cor- 
related with V (Figure 5) but this was probably an instrumental effect. 
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and the dispersion equation is 


0)2 _ k 2 v^ 2 cos 2 0 - 0. (4.20) 

This is the only mode that can propagate in an incompressible hydromagnetic 

fluid. 

Using an iteration technique, Unti and Neugebauer showed that the observations 
were consistent with the condition n = constant (see Figure 5). 

Unti and Neugebauer (1968) searched for sinusoidal Alfven waves in the solar 
wind using magnetic field and plasma data from Mariner II. One result is given 
in Figure 4, This shows B (smoothed by a 19 point running average) in a co- 
ordinate system for which the minimum field fluctuation was in the xy plane and 
the average field along z, the bulk speed of the plasma is also shown. One sees 
a nearly sinusoidal, plane polarized wave with 8 b strongly correlated with V 
and S B normal to B. Using the measured B and n to compute Sv they found good 
agreement with the measured bulk speed perturbations (Figure 5), suggesting 
that (4.17) is satisfied. The initial density measurements, computed on the 
assumption of a strictly radial bulk velocity, did not satisfy the condition n = 
constant, necessary for an Alfven wave, (in fact the density was strongly cor- 
related with V (Figure 5) but this was probably an instrumental effect. 
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Belcher et al. (1969) reported that correlation coefficients >.8 in absolute 
value between the radial component of the interplanetary magnetic field and solar 
wind speed were found at least 30% of the time in the 5 months of data from 
Mariner 5. They argued that these are due to large amplitude, aperiodic Alfven 
waves because: (1) the correlation is consistent with (4.17), (2) the density was 
not correlated with the field changes, (3) the value of B/V a computed from plots 

i 

of b versus V agrees with that given by (4.17) with £ = 1 (see Figure 6a) , and (4) the 
fluctuations SB tended to be normal to B (see Figure 6b). They also found that 
essentially all of the waves were propagating away from the sun, irrespective 
of the direction of the interplanetary magnetic field. This is evidence for the 
existence of Alfven waves in the solar wind and suggests the relatively common 
occurrence of such waves on a scale of several hours or more. 

Coleman (1966) had previously noted a high value for the coherence he obtained 
from the cross spectrum for V and B , and he noted that the phase difference at 
a particular frequency obtained from the cross spectra for V ^nd b was ~ 180° 
when B was away from the sun and ^0° when B was toward the sun. These 
evidently correspond to the results seen by Belcher et al. (1969), expressed 
in a different form. Coleman concluded that fast waves or Alfven waves (or 
both) were present, but he could not distinguish between the two. A more 
detailed account of this work also does not give positive evidence for Alfven 
waves (Coleman, 1967). 
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2. Magneto Acoustic Wave (G ^ 0, M = 0; Vj| = 0). While the Alfven 


waves propagate along B because of the tension in the field lines, the restoring 
force for these magneto acoustic waves arises from the compression of p and 
B , and the wave propagates normal to B (see Figure 2b). Taking M= (B • V ) 

B = 0 (i.e., assuming the field lines do not bend) and B - (B + b z ) £ , the linear- 
ized form of (3.12) becomes 


P 


B v 
B~t 


i 




(4.21) 


for a wave propagating along x. The elements of the wave are shown in Figure 
7. The frozen-in condition (3.13) gives 


3 V„ . 


Bb z _ ^ 

B t B x 


(4.22) 


Operating on (4.22) by B /B t and using the continuity equation gives 


B 2 b 
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Bt 2 


- ( c l + V 2 ) 
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where 


c 2 ~ ffX 


H P- L 


and c 2 ,' 


3P, 




(4.23) 


(4.24) 


Thus, the wave propagates along x with the speed 


V 2 = c 2 ^ V 2 
v MA c l T V A’ 


(4.25) 


The dispersion equation is 
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c 2 - k 2 (c 2 + V 2 ) = 0 


(4.26) 


Burlaga (1968) has presented evidence for a magneto acoustic wave in the 
Pioneer 6 magnetometer data of Ness. This is shown in Figure 8 together with 
Pioneer 6 plasma data from the MIT probe. One sees sinusoidal oscillations in 
I B |, but none in B/B, as expected for a magneto acoustic wave. The plasma 
density and temperature should oscillate in phase with B . The density is given 
in Figure 8 to the nearest integer, reflecting the accuracy; this "is not adequate 
to show the small oscillations that are expected. The temperature measurements 
suggest a maximum at the time of maximum compression, but they are not ac- 
curate enough to establish this with certainty. Note that to resolve these waves 
with a plasma probe it would be desirable to measure at least 2 spectra per 
minute, and since A n/n % AB/B = .1, the density should be measured with a 
relative uncertainty appreciably less than 10%. 


8. lT Fast Waves" (G ^ 0, T ^ Q; v n = 0) . In this case both tensile and 
compressive stresses are active, but there are no displacements along B . Let 
x be along § v and take k and S B in the x - z plane, (see Figure 9). The pro-- 
cedure of the preceding 2 sections gives the following equation 


3 2 b 
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3 t 2 


(Cf + v 2 ) 


B 2 b 
B x 2 
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B 2 b 
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B z 2 


- ,0 


(4.-27) 
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and a similar equation for v x . The dispersion equation is 


- (cl + V 2 ) k 2 _V' A 2 k 2 = 0. 


(4.28) 


Clearly, the mode is related to the Alfven wave and the magneto acoustic wave 
discussed above. 


If the compression is adiabatic with an adiabatic exponent y = 2, then = 
A i p 2 c 2 = 2 2 p and (4.28) can be written 


co 2 - a 2 k 2 - b 2 k 2 = 0 

Z X 


where 


and 


a 2 



2 P, )/P = V 2 + 


(P L ~ Pi]) 


b 2 B 





2 Pi 


P 


(4.29) 


(4.30) 


(4.31) 


Equation (4.29) is precisely the equation derived by Whang. (1971) for a guiding 
center plasma. 


Let the angle between k and B be 8 . Then k~ = k 2 sin2£and k 2 = k 2 cos e, 
and the dispersion equation (4.28) has the form 


<o 2 _ k 2 [(c 2 + V 2 ) - (c 2 + (1 - a) cos 2 0=0 ( 4 * 32 ) 
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This reduces to the dispersion equations derived above for Alfven wavbs and 
magneto sonic waves when <9=0 and n /2 (equations (4.20) and (4.26)), 
respectively. 

Evidence for this type of wave was given by Whang and Ness (1-970). They point 
out that the moon will generate such waves when the solar wind moves past it, 
as a result of the different magnetic permeability of. the- moon and the solar 
wind. The waves will form an elliptical' Mach cone with the long, axis perpen- 
dicular to the average 1 field direction and the short axis along B;, as a result of 
the asymmetry of the- hodograph corresponding to (4.32). Crossings! of the 
cone were identified in the magnetic field- data from Explorer 35, as described 
by Whang and Ness (1970), and they obtained the Mach cone cross-section shown 
in Figure 10'. It is an ellipse whose ratio of major and minor axes is in agree- 
ment with the fast wave theory given above. Michel (1967) had' earlier predicted 
an elliptical Mach cone on the basis of a qualitative extension and application of 
ordinary, isotropic gas dynamics. 

4. Coupling Between the Fast Wave and the Alfven wave. We have 
shown that a "fast" wave can propagate with k, 8 B and § v in the x-z. plane. 

This becomes an Alfven wave when k is along B and a magneto acoustic wave 
when k is normal to B . In addition, an Alfven wave can propagate along k , with 
8 v and SB along y , normal to the x-y plane. Thus in general both types of waves 
can propogate simultaneously. Their hodographs are shown in Figure 11. 
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Iti the linear approximation, the fast wave and the Alfven wave are uncoupled 
and propagate independently.; The fast mode causes oscillations in the intensity 
of B Both the fast mode and the Alfven mode cause fluctuations in the direction 
of B ; in general, these polarization of these oscillators will be elliptical. 

Burlaga found a wave in the Pioneer 6 data which has the characteristics just 
described (Figure 12). Again, however, the observations were not complete 
enough or accurate enough to verify the theory just presented. 

The fast wave and Alfven wave in Figure 11 are uncoupled on the linear approxi- 
mation. However, if the amplitudes of the oscillations are large, one must 
consider the effect of higher order terms. If the term of second order in b 2 
is not discarded in the compressive term of the force equation (3.12), one gets 
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3t 2 


- (cf + V 


2 ) 

A 1 


h 2 b 

x 
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3 2 (b 2 ) 

3 x 2 


(4.33) 


which reduces to (4.27) when the term in b 2 is small. Thus the fast wave, which 
was confined to the x~z plane, can couple with the Alfven wave that has perturba- 
tions along y , through the non-linear term on the RHS of (4.33). Kawashima 
(1969) first noted this (but there are some minor errors in his equation cor- 
responding to (4.33)). He suggested that magneto acoustic waves are produced 
by non-linear Alfven waves in this way in the solar wind. 
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The power observed in the compressive waves is always lower than that in 
Alfven' waves for frequencies in the range 10" 5 to 10 “ 2 Hz -(e.g. Coleman, 1966) 
which is 'consistent with Kawashima's theory. Clearly, the coupling should be 
larger the larger the amplitude of the Alfven waves. Thus, the theory implies 
that the proportion of compressive waves to be large when the amplitude of the 
fluctuations is large. Kawashima used Ness' IMP 2 observations to show that 
this is the case {Figure IS). 


5. Firehose Instability. The dispersion equation for Alfven waves 
(4.14) shows that waves grow with time {co is imaginary) when < 0., i.e., 
when 



B 2 /(8 t t) 


< 2 


(4.34) 


Parker ,(1957) was the first to suggest that this instability might occur in the 
solar wind. He derived an equation which is the same as 4.34 with i = 1, i.e., 
he considered protons only. 


Kennell and Scarf (1968) pointed out that electrons should also be considered, 
and they obtained (4.34) with i = 2. Kennel and Scarf derived their result from 
the Boltzman equation. They wrote (4.34) in the form 


2 - A, A* - J3_ A_ < 0 


(4.35) 
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where 


and 


/3± S 


p \]± 

B 2 / (8 it) 



(4.36) 


(4.37) 


The growth rate for the firehose instability is characterized by y where 



(4.38) 


Eviatar and Shultz (1970) found that the growth rate is so large in the solar wind 
that one would not expect to find 


(B 2 /8 t r)' 1 ^ (Pin - P u ) >2.1. 

i 

in the solar wind, i.e., the firehose instability, would limit the thermal anisotropy. 
This raises the question, Does the anisotropy ever get large enough in the solar 
wind for the firehose instability to grow, or does some other mechanism destroy 
the anisotropy first? 

There is no unambiguous evidence for the existence of the firehose instability in 
the solar wind atlAU. Scarf et. al. (1967) using average values for the solar 
wind parameters, concluded that (4.34) for protons is not satisfied, so the fire- 
hose instability is not expected to develop. Kennel and Scarf (1968) argued that 
protons alone cannot cause instability. But they argue that since the field is 
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disordered the electrons must be sufficiently hot and anisotropic that (4.34) is 
satisfied (assuming that the firehose instability is the cause of the fluctuations). 


The most recent report of measurements by the Los Alamos group (Hundhausen 

et..al., 1970) gives T + /T + = 1.9, T a +„ - 7.2 x 10* °K, and n 

7.7 cm" 3 . The electron temperature seems to be essentially constant, T ~ 

(1.5 ± .5) x 10 s °K, as discussed in Section II, and the best measurements of 
the electron anisotropy give .1 ^ (T|~ - Tj^)/T" < .2 . The average field is 
typically ^ 6 y. These numbers give 2 -/3 + x A + - j3_ x A_ = .8, so again one 
concludes that the firehose is not normally operative unless the electron ani- 
sotropy is larger than the current measurements indicate. One factor has still 
been left out - the a- particles; but their anisotropy would have to be very large 
to have a significant effect, since their density is relatively low. 


Eviatar and Shultz (1970) used Vela data to compute - /3^ for a single 5- 
hour period. The results (Figure 14) show that (4.34) was not satisfied for that 
period. 


Burlaga et. al. (1969) considered the possibility that the instability might develop 
in some regions but not in others. They identified very quiet and very disturbed 
one hour intervals in the Explorer 34 ma gnetic field data (Figure 15) and com- 
puted /3 + for each of the intervals using the plasma and magnetic field data from 
that spacecraft. These intervals tended to be isolated, indicating that the dis- 
turbances were highly localized. They found that /3 + tended to be high for the 
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disturbed intervals, but low for the quiet intervals, which is, of course, what 
(4.34) would imply. Assuming that .1 4 (T“ - )/T | |~ < .2, they found that A + 

+ KA“ < 1.5, where K = T"/T + , so the instability condition (4.34), which can be 
written /3 + (A + + KA~) < 2 is satisfied only if > 1.3. Figure 15 shows that in fact 
the largest percentage of very disturbed intervals did occur when /3 >1.3. Thus, 
there is evidence that necessary Conditions for the firehouse instability do occur 
occasionally. However, there were .very disturbed conditions when /3 + was less 
than 1; in fact the most probable /3 + for the disturbed intervals was .7. The 
most probable value of f3 T for the Explorer 34 data (/3* p = .3) is appreciably 
less than that found by Neugebauer and Snyder (1967). The differences could 
be due to a systematic error in the Explorer 34 density measurements; but the 
absolute accuracy of density determination was estimated to be ±20% (Ogilvie 
et al. (1967), and this is consistent with other observations (Burlaga and 
Ogilvie, 1970; Gilbert, Private Communications.) 

Summarizing, the firehose instability probably does occasionally occur in iso- 
lated intervals where /S + is large and possibly for most "very disturbed intervals", 
but it does not seem to be the primary cause of the fluctuations that are 
typically observed. Clearly, more measurements are needed. 

Belcher and Davis (1971) reported a poor correlation between /3 + and cr 
the square root of the 3- hour average of the 168.75 sec minute total variance 
in field components. They considered this to be in conflict with the results of 
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Burlaga et al. (1969). There is no conflict, however. cr gi i S not a good measure 
of the type of high frequency fluctuations in the magnitude and direction of 8 
studied by Burlaga et al.. 

D. V|j ^ 0, v± 0 (General MHD Dispersion Equation), We have been dis- 
cussing various special cases of low frequency, long wavelength waves. In gen- 
eral, all of the modes might exist and might be coupled. From isotropic mhd 
wave theory (e.g., Jeffrey and Taniuti (1964)) it is well known that 3 modes would 
then propagate: (1) the fast magneto acoustic mode, which reduces to the "fast” 
magneto acoustic wave discussed above for propagation normal to B , (2) The 
Alfven mode, and (3) the slow magneto acoustic mode, which reduces to the ion 
acoustic mode for propagation along B. The dispersion equation for an anisotropi 
medium is very complicated (It is given in Lust, 1959), but gives the same quali- 
tative features as that for an isotropic medium. The surfaces of normal speeds 
(hodographs) are shown in Figure 16. The fast magneto acoustic mode is the 
one most analagous to the sound wave in gas dynamics, in the sense that it 
propagates in all directions. Here, however, the propagation is anisotropic. 

The Alfven mode cannot propagate normal to B but it does propagate in all other 
directions, the speed being greatest and equal to the Alfven speed for propagation 
along B . The slow magneto acoustic wave has characteristics of both the Alfven 
wave and the sound wave; it does not propagate normal to B, but it does propa- 
gate in all other directions, the speed being greatest for propagation along B . 
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Figure 16 shows that 3 ho dographs must be distinguished, corresponding to 
V' < c.. , Va' < c.. and V ' = e . Although a distribution of V' / c has not been 

a II SI A H A [| 

published, it is clear from the measure parameters that V^/c'is typically very 
close- to 1. This is an important point which has not been discussed enough in 
the literature. In principal, one can excite all 3 modes simultaneously by an 

V 

instantaneous point disturbance in the solar wind, assuming that the ion acoustic 
wave propagates and is coupled to other modes . It is of interest to ask what is 
the form of the 3 corresponding wave fronts. For the three cases correspond- 
ing to the hodographs in Figure 16 the fronts at any instant are as shown in 
Figure '17. At a later time the fast wave front is larger. The slow wave front 
also enlarges with time, and it propagates along 3 • 

As discussed earlier, the moon can act as a source of waves which propagate 
away from the moon while being convected by the solar wind. If the above theory 
applies, these waves will form one Mach cone corresponding to the fast mode, 
which is analogous to the aerodynamic Mach cone and equivalent to that pre- 
dicted by Whang's model, and two Mach cones corresponding to the slow modes, 
when V [ c.. . (When V' = c , there is just one Mach cone.) These can be 

constructed from the Friedrichs diagram as illustrated in Figure 18. Under 
typical solar wind conditions near 1 AU it would be very difficult to observe the 
slow Mach cones, and in fact they would not always exist. 

There is only indirect evidence for the existence of the slow mode on the solar 
wind. This is the observation of slow shocks, discussed in Section V, 
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V. Pis continuities 

The theory of discontinuities in a steady state, isotropic hydromagnetic fluid is 
described in many texts (e.g. Jeffrey and Taniuti, 1964) and reviews (Colburn 
and Sonett, 1966) but the theory for an anisotropic fluid has been developed only 
recently in response to solar wind measurements. Lynn (1967) has presented the 
theory for a single species plasma in the CGL approximation. A similar theory 
for shocks was given by Abraham-Shrauner (1967). A less restrictive theory for 
shocks was given by Chao (1970) and a similar theory including other types of 
discontinuities was published .by Hudson (1970). Neubauer (1970) has presented a 
theory of shocks in an anisotropic plasma which parallels that of Chao and Hudson. 
Our treatment will most resemble those of Chao and Hudson. 

Let the discontinuity be a plane surface which moves with speed U. along its 
normal n relative to a fixed frame (the sun, say). Let the subscripts 1 and 2 
denote the regions away from the sun and toward the sun, respectively. Let the 
subscript t denote components along some direction parallel to the discontinuity 
surface. 

Assume a steady state and work in a frame moving with the discontinuity surface. 

Then, applying Gauss’ theorem to (3.1) shows that the mass flux normal to the 
surface, G, is continuous across the surface, 

. G s Pt 0 In - U) = P 2 ( v 2n - u )- t 5 ' 1 ) 
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Since V- B = 0 


B *B. . 

n In 2n 


(5.2) 


The frozen- field condition (3.13) implies that 


B n <* U “ V 2t^ = B lt ^ v in ~ B 2t ( V 2 n " U ) 


(5.3) 


The momentum flux equation (3.12) gives 2 relations, one for the change in. the 
normal momentum flux 
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another for the change in the transverse momentum flux 
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Using (5.3), this can be written 
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which with (5.1) gives 
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where 


and 



B 2 

n 

AlT ~P\ 





B 2 

n 

4 TT P 2 


^2 


are the Alfven speeds along the shock normal. 


(5.8) 


Finally, the energy equation (3.22) gives the following relation for the change in 
the kinetic energy flux normal to the surface. 



e +E PkL 
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“ B n B t ' v t ^/( 4 ^ 


where e is given by (3.18). 

Equations (5.1), (5.1), (5.3), (5.4), (5.5) and (5.9) are the basic jump equations. 

There are 8 equations for the 9 variables b , v , p, p and p . 

~ ~ l II 

Using (5.1) and (5.2), equations (5.3), (5.4), (5.5) and (5.9) can be written as 
follows 


G 




(5.10) 
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which, with (5.1) gives 
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(5.14) 


Two general classes of discontinuities can be distinguished: 


(1) Stationary discontinuities (G - 0) 

(2) Shocks (G ^ 0) 

B. Stationary Discontinuities. Since by definition G = 0 for a stationary 
discontinuity, v nj = v n2 = U, which implies that the surface does not move rela- 
tive to the plasma. There are two types of such discontinuities, corresponding 
to B n =0 and B n / 0. These are called tangential discontinuities and contact 
discontinuities, respectively. 

1. Tangential Discontinuities (G = 0, B n - 0). There are several 

necessary conditions for a tangential discontinuity: 
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(a) Since G = 0, it is non-propagating, i.e., the speed of the dis- ' 
continuity surface relative to the solar wind in the direction of its normal is 

* r ■> 

zero. 


(b) B n = 0 implies that ,B 1 and B 2 are parallel to the discontinuity 
surface {see Figure 19); hence the normal is 


n = B 1 x B 2 /| B 1 x B,_ 


(5.15) 


(c) Since there is no change on V n , v, and v; 2 are parallel to the 


surface; hence 


n = v x x v 2 /| Vl x 


(5.16) 


(d) The equation for the normal momentum flux (5,12) gives 




= 0 , 


(5.17) 


which says that the pressure perpendicular to the surface is the same on both 


sides of the surface. Since - n k * 1^ k ’ this ma y wr hten 




= 0 . 


(5.18) 


The equations for' the transverse momentum flux and the energy flux, (5.13) and 
(5.14) are identically satisfied. 
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To positively identify a tangential discontinuity, one must show that all of the 
above conditions are satisfied, which requires simultaneous magnetic field data 
( B (t)) and plasma data(n, 1) , and v for protons, a’s and electrons) and the ar- 
rival times at 4 spacecraft (to get the speed and oiientation) . Such a definitive 
measurement is still well beyond the scope of the observations, but there are 
several somewhat less complete results that are consistent with these conditions 
and leave little doubt about the existence of tangential discontinuities in the solar 
wind. 


Figure 20, from Burlaga (1968), shows a region with at least 3 discontinuities 
which were observed by the spacecraft Pioneer 6. The discontinuity surfaces 
are moving past the spacecraft, nearly radially away from the sun with a speed 
' Xj 400 km/sec. They are seen as abrupt decreases in the magnetic field intensity 
(from Ness’ magnetometer experiment) and simultaneous increases in the proton 
density (from the MIT plasma detector). The change in B occurs in£ 30 sec, 
implying a thickness -$12,000 km « 10 4 AU, assuming V - 460 km/sec, which is 
indeed discontinuous on a scale of 1 AU. If the region in Figure 20 was iso- 
thermal and if n a = an p there, then the total pressure is constant throughout the 
region if and only if n ~ B 2 , since 
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+ B 2 /(877) 
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Figure 20 shows that in fact n is proportional to B 2 for the observations be- 
tween 0630 and 0730, Thus the data are consistent with a constant pressure 
across each of the discontinuities, in agreement with (5.17). 

Condition' (5.15) can be tested only for discontinuities across which there is a 
significant change in the direction of J5 . Burlaga (1969a) has introduced the 
terms "directional discontinuity" for discontinuities across which the field 
direction changes by> 30° in < 30 sec. Figure 21 (from Ness et. al., 1966) shows 
2 directional discontinuities, one seen at Pioneer 6, the other seen at Explorer 
33. 


Computing the orientation of the discontinuity surface from (5.15) and using the 
measured -solar wind speed 400 km/sec, Burlaga and Ness (1969) computed that 
the discontinuity should have arrived at IMP 3 57.5 min after it was seen at 
Pioneer 6. The observed transit time was 53 min. Thus, the observations are 
consistent with the assumptions that B was parallel to the surface and that the 
surface was convected with the solar wind speed, as required for a tangential 
discontinuity. 

Burlaga and Ness (1969) used simultaneous magnetic field measurements from 
3 spacecraft (Explorers 33, 34, and 35) to identify 6 directional discontinuities 
whose orientations determined from measured time delays and solar wind 
speeds, are consistent with those computed from (5.15). These results are 
shown in Figure 22. Consider panel A in this figure. A discontinuity was 
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observed first by the spacecraft at the point marked 1. Using the magnetic field 
data from that spacecraft, the surface normal was computed assuming (5.15) 
and the intersection of that surface with the ecliptic plane was computed and is 
shown by the line passing through point 1. A few minutes later, the same dis- 
continuity was seen at another spacecraft at the point marked 2. The surface 
intersection was similarly computed using data - from the second spacecraft at 
position 2. Minutes later it arrived at the spacecraft at point 3. Each 
of the panels refers to a different discontinuity surface and has a 
similar interpretation. If the surfaces are tangential discontinuities and if 
they are plane, the lines through points 1, 2, and 3 on each panel should be 
parallel. This is the case in A, D, E, and F, but in B and C the surfaces seem 
to be appreciably curved on a scale of 50 R_. Note that the order of arrival of • 
the surfaces at the different S/ C depends on the surface orientation computed 
from (5.15). In every case, the arrival order is consistent with the assumption 
that the surfaces are tangential discontinuities with orientations given by (5.15). 
In this connection, D, E, and F are particularly noteworthy. Burlaga and Ness 
(1969) have shown that the actual surface orientations determined from the delay 
times and the measured solar wind speeds, are in very good agreement with the 
theoretical orientations shown in Figure 22. They also showed that (5.15) is 
satisfied for these discontinuities. ■ 

Equation 5.10 shows that^ , and v 2t are arbitrary for a tangential discontinuity, 
since G = 0 and B n = 0. (This does not consider the effect of instabilities which 
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might arise). In particular, since v Jn = v 2n > jv u j ^ | v 2t | implies that one may 
have- ‘v , ^ {v j.-. Such changes in the bulk speed do occur at directional dis- 
continuities (Burlaga, 1968). Burlaga (1969) and Siscoe et al. (1969) have found 
discontinuities across which the bulk speed changed by 60 km/sec or more in 
less than ^3 min. 

Siscoe et al. (1969) pointed out that in addition to conditions (5.15) to (5.18), a 
tangential discontinuity must satisfy the condition 

A v . 

tan 0 = s tan iir (5.19) 

A v k 

where the quantities 0 , Av t and Av k are as defined in Figure 33 a. This simply 
says that a unit volume element initially adjacent to a tangential discontinuity 
surface always moves parallel to this surface as the surface moves away from the 
sun. If a discontinuity is tangential, 0 can be computed from (5.15) using 
measurements of B, and it can be computed independently from (5.19) using 
measurements of v . Siscoe et al. (1969) did this for 12 discontinuities in the 
Pioneer 6 data and found the results shown in Figure 23b. There is considerable 
scatter because of the relatively large errors in measuring the flow direction; 
in fact the error bars (not shown) are so large that only 1 point deviates signifi- 
cantly from the line 0 = 0 . Siscoe et al. concluded that the discontinuities were 
tangential. 
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We have been discussing discontinuities as if they have zero thickness, i.e. true 
discontinuities in the mathematical sense. Physically, of course, the surface 
has some structure. This has lead to some confusion. Figure 24 shows the 
"thickness" distribution of "discontinuities" in the magnetic field, studied by 
Siscoe et al. (1968). The "thickness" is here measured by the time for the sur- 
face to pass the spacecraft, which is on the order of the actual thickness divided 
by the solar wind speed. The shape of the distribution is very important, for it - 
shows that these discontinuities are indeed a distinct class of thin structures .and 
not simply the tail end of the distribution of fluctuations seen in space. Since the 
directional discontinuities studied by Burlaga (1969a) has the same statistical 
properties as those studied by Siscoe et al. (1968), they probably have the same 
kind of thickness distribution. The discontinuities studied recently by Smith et al., 
(1970) and Turner and Siscoe (1971) have T> 300 sec, which puts them far to 
the right of Figure 24 and may indicate that they are a different class of objects. 

It is not clear whether the structure of tangential discontinuity surfaces in the 
solar wind can be analyzed within the framework of magnetohydrodynanrics. In 
any case, it has not yet been done. One expects the fluid theory to be applicable 
on scales appreciably greater than the proton gyro radius, corresponding to T» 

1 sec, which seems to be the case for most of the surfaces in Figure 24. Two 
types of structures have been identified in the magnetic field: (1) a laminar cur- 
rent sheet, in which 8 rotates more or less uniformly in the plane of the dis- 
continuity surface as shown in Figure 25 (Siscoe et al. (1968), Burlaga and Ness 
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(1969)) and (2) a "turbulent" structure, as shown in Figure 26 from Burlaga 
(1969b). The study of such structures is still in an early stage. 

Northrop and Birmingham (1970) studied the stability of discontinuity surfaces 
of finite thickness. Using an MHD energy principal due to Grad and Rebhan 
(1969), they gave an elegant proof that such surfaces are stable if v lt = v 2fc , 

In addition to the structure of the current sheets associated with tangential dis- 
continuities, illustrated in Figure 26, there is another type of structure which 
Burlaga and Ness (1968) called D- sheets. Two particularly noteworthy examples 
are shown in Figure 27, from Burlaga (1968). These are characterized hy a dip 
in the intensity of B over an interval of a few minutes associated with a discon- 
tinuity in the direction of B. Several other examples are shown in Figure 28 
from Burlaga and Ness (1968). The physical significance is not understood. If 
the dip were caused by the annihilation of the antiparallel components of B l and 
B 2 , then one expects the minimum field to be related to the angle co between Bj 
andB 2 , by the relation B m i J B = cos &>/2, if B x = B^. Figure 29 (bottom) 
from Burlaga (1968) shows that this is the case for the D-sheets discontinuities 
in Figure 28 (top). It is not certain that the discontinuities with = B 2 are 
tangential. Figure 27 shows one D sheet associated with a large change in B 
and n, which is almost certainly a tangential discontinuity; here too, is 
consistent with the annihilation hypothesis. 
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2. Contact discontinuity. In this case G = 0 but B ^ 0. It 

n 

is convected, but (5.10) shows that in this ease one must have v =v 2 t» i,e -> 
there can be no relative motions along the discontinuity surface. This is be- 
cause 8 is frozen to the plasma; a shear when B n £ 0 is impossible without 
breaking lines of force. In an isotropic plasma, conservation of transverse 
momentum (5.13) requires that B lt - B 2t . So there is no change in J3 across a 
contact discontinuity in an isotropic plasma— no change in the magnitude, no 
change in direction (see Figure 29). In an anisotropic plasma, however, one can 
haveB lt ^B 2t if^j ^ (see (5.13)). Thus, in an anisotropic plasma one can 
have a change in the magnitude and direction of b across a contact discontinuity, 
but here as in the isotropic plasma v lt = v 2t . 

The pressure balance (normal momentum flux) condition (5.12) is 

b2 B 2 1 2 

£ p kl + 87 - 57 I>k» - Pkpj = 0 (B.20) 

1 

This has a simple geometrical interpretation. Let 6 be the angle between B and 
n. Then the pressure normal to the surface is 

p n = p !i n + p i n = P|| - cos2 9 + (pi + 1^) sin2 9 > < 5 ‘ 21 > 

where 

B 

cos 6 — — 

B, 
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Setting p nl = p n2 gives (5.20). Contact discontinuities have not been observed 
in the solar wind. It is widely believed that they will not be observed because 
they would rapidly broaden into a smooth transition. 

C. Shocks 

1. Introduction. A discontinuity is called a shock if G £ 0. In aero- 
dynamics there is basically just one land of shock, and it correspond to the one 
type of wave in an ordinary gas, the sound wave. Gas enters an aerodynamic 
shock at a speed greater than the sound speed and leaves at a speed.less than 
the sound speed. The energy is lost from the streaming goes into heating the 
gas. In magneto-gas-dynamics there are 3 basic wave modes, as discussed in 
Section IV and there is a corresponding variety of shocks. More complicated 
energy transfers are possible because of the presence of the magnetic sink. 

We shall consider first the simple cases of shocks moving parallel to B or 
perpendicular to B . Then we consider shocks more generally, shocks moving 
in arbitrary directions with respect to B . Some special types of shocks will be 
discussed at the end of this section, together with a classification scheme. A 
detailed analytical theory of hydromagnetic shocks can be found in. Jeffrey and 
Taniuti (1964). 

In principle, a shock in the solar wind can move either away from the sun or 
toward it* Shocks moving away from the sun are called forward shocks, while 
shocks moving toward the sun are called reverse shocks. 
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Physically, there is no difference between a forward shock and a reverse shock, 
but there is probably a difference in the origin of the forward and reverse shocks 
that are observed in the solar wind. The origin of the observed reverse shocks, 
described below, is still a mystery, 

2. Parallel Shocks. A parallel shock is one whose normal is parallel 
to B . (See Figure 30). Two waves can propagate along B, the acoustic wave 
and the Alfven wave. So immediately the question arises, "with respect to which 
of these modes is the shock supersonic?" The answer is somewhat complicated. 

We shall just present the results of Jeffrey and Taniuti (1964, p. 246). 

If the sound speed is greater than the Alfven speed .ahead of the shock, then the 
sound speed is dominant and one has an ordinary aerodjmamic shock with gas 
entering faster than the sound speed e N and leaving at a speed less than c ([ . The 
flow is super- Alfvenic on both sides of the shock and the sound speed exceeds the 
Alfven speed on both sides. 

If the Alfven speed is larger than the sound speed ahead of the shock, there are 
3 possibilities. The obvious one, with gas entering super-Alfvenically and 
leaving sub-Alfvenically probably does not occur in nature (Jeffrey and Taniuti 
1964). Another possibility, with gas entering supersonically (and super- 
Alfvenically) and leaving subsonically (but still super-Alfvenically), is allowed 
and might be observed. Finally, gas might enter supersonically (but sub alfvenic ally 
and leave at a speed which is either sub Alfvenic or subsonic. 
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For parallel shocks, (5.10) implies that v = v . There is no requirement 
that v i , v 2 > n be coplamar. The magnetic field drops out of the Rankine 
Hugoniot equations, which become identical to those for aerodynamic shocks, 
except that they include both P ± and p ,, in the expression for internal energy, e . 
These equations provide no insight concerning the possible change in anisotropy 
(distribution of internal energy) across a parallel shock. 

There are no observations of parallel shocks in the solar wind. Such observa- 
tions would be of special’ interest. 

3. Perpendicular Shock. A shock is called a perpendicular shock if 
the shock normal is perpendicular to 8 . In this case B n = 0 and the magnetic 
field is parallel to the shock surface on both sides of the shock, as shown in 
Figure 31. 


The only wave that propagates perpendicular to B is the magneto acoustic wave, 
whose speed is V M 2 = V 2 + Cj 2 . The flow enters a perpendicular shock at a speed 
greater than V M and leaves at a speed less than V M , Thus, this shock propagates 
supermagnetoacoustically. 


For a perpendicular shock, the mass flux and frozen field equations (5.11) give 


Py Pn 


(5.22) 


which states that (a) the magnetic field direction does not change across the shock, 
and (b), its intensity changes in the same ratio as the density, B /JB^ = p 1 /p 2 . 
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Equation 5.13 shows that there is no change in the transverse momentum flux 
across a shock when B n = 0, so v lt - v 2t . This implies that v 2 - v 1 = 

(v 0n - v in )h,so that h is parallel to Vl - v 2 for a perpendicular shock. 

The normal momentum flux equation (5.4), gives 

P (V n - U)2 

i.e. the change in the normal momentum flux is balanced by a change in the 
pressure. This is the same as for an ordinary gas dynamic shock, except that 
here the magnetic field pressure adds to the particle pressure. 

Equation (5.9) shows that the change in the kinetic energy flux normal to the 
surface is 



which is likewise the familiar gas dynamic result with the addition of the mag- 
netic pressure to the thermal pressure. 

Note that the parallel pressure does not appear in the Rankine-Hugoniot equations. 

4. Reverse Perpendicular Shock. As we have said, a reverse shock is 
me which propagates toward the sun relative to a frame moving with the plasma. 
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Since the solar wind speed is usually higher than the shock speeds, the shock 
can actually move away from the sun while propagating toward the sun. In this 
case, a spacecraft sees the "back," high entropy, side first and then the "front." 
This means that the density, temperature, and magnetic field intensity should be 
seen to decrease simultaneously as the shock is carried past the spacecraft. 

As always, the speed behind the shock is lower than that ahead of the shock, so 
in the case of a reverse shock ( v nl ~ U ) ( v n2 - JU )» since the subscript 1 refers 
to the side away from the sun. Thus, > v nl and an observer at a fixed 
spacecraft would see the bulk speed increase. 

While- a reverse shock looks different in the data in the sense just, discussed, it 
satisfies the same equations as the forward shock (taking proper account of the 
sign of U, of course). 

A perpendicular shock was found in the solar wind by Ogilive and Burlaga (1969) 
using Explorer 34 magnetic field and plasma data. It passed the spacecraft at 
1732 UT on Aug. 29, 1967. As shown in Table III, the bulk speed, density, 
temperature and magnetic field intensity all increased, simultaneously, but there 

t 

was no change on the direction of the magnetic field. The shock surface moved 

past 3 spacecraft, each carrying one of Ness 1 magnetometers. The transit times 

* 

were measured, and using these the orientation of the shock surface was com- 
puted with the result shown in Figure 32. The shock normal determined in this 
way is n (-.93, -.33, 1.5) in solar ecliptic coordinates i.e., & n = 9°, 0 n = 200°. 
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Table III 


1732 UT, Aug. 29, 1967 Shock 



V (km/sec) 

n(eru” 3 ) 

T(°K) 

B(7) 

e . 

(p 

Before 

418 

2.6 

6.5 X 10 4 

5.5 

51° 

295° 

After 

452 

3.7 

12 x 10 4 

7.2 

51° 



296° 
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TMs is essentially normal to the magnetic field, B 1 and B 2 , as required for a 
perpendicular shock. 

From the multi spacecraft measurements, the speed of. the, surface was found to 
be U = (436 ± 110) km/se.c. This compares favorably with . the value obtained 
from the mass flux conservation condition (5.1), viz U = 496 km/sec. 

Equation (-5.22), which -is a consequence of the frozen field condition, is 
satisfied -since, (a)’ the field direction did not change, and (b),the measured ratio 
Bi/ n 1 = 2.0 equals the ratio B 2 /n 2 = 1.9 within the experimental. errors. 

The change in the* momentum' flux, p ,(V n ~ U) j , was 1.4 x 10“ 10 dynes/cmd , 
which is equal to the change in pressure 



+ 


52 

8 7 7 


2 

= 1.5 x 10“ 10 dynes/cm 2 * * 


1 


(The pressure .here was computed assuming that' T j = T = 1.5 x 10 5 °K and 
that the contribution of the a's. is negligible.) Thus (5.24) is satisfied. 


Finally, the change in the kinetic energy flux given by the LHS of (5.25) was 
consistent with the observed change in the RHS of (5.25) within the errors. 


Burlaga (1970) identified a reverse perpendicular shock in the Explorer 34 
plasma and magnetic field data of Ogilvie and Ness. The observations are 
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Table IV 


Sept. 28, 1967 



V (km/sec) 

n (cm. -3 ) 

T (°K) 

B(r) 

0 

CP 

Before 

540 

10 

4.5 x 10 5 

12 

© 

o 

130° 

After 

585 

5.5 

2.6 x 10 5 

8. 

0 

CO 

125° 

















shown* in 'Figure 33 where^he- decrease: -in ; n„ T, and B,. and .the.fnerease in V 
required for a reverse perpendicular shock are clearly seen.. Table-- IV shows 
the -magnitude of these changes- together with the -magnetic 'field- dire cti ons . The . ' 
magnetic field direction did not change, within the measurement 'uncertainties ' 

(= ± 5°), as required for a perpendicular shock (Equation (5.22)). 

The- shock. surface moved past 3 spacecraft,, each carrying , one; of ',-Nes s' ! mag- 
netometers, just as the forward perpendicular shock. discussed-in the previous 
section. The shock normal was foundfo be n = (-.70, - . 69,, ‘.18,), i.e., 9 = 10% 

n 

4> n - 225°. This is -essentially normal to the magnetic field-fcos -1 (h • B.) = 87°), 

‘as required for a -perpendicular shock. 

Using the multispacecraft .'magnetometer-measurements-, of ISfess ,and : the- conser- 
vation of mass flux equation (5.1), the -shock surface was found' -to -be 
moving at a speed U = 2-70 km/sec. in .the* sunward direction. Since. this 
speed is less than the component of the solar wind speed. ,along*the 'normal, the 
shock was carried away from the sun- despite its motion-. toward the ‘ sun relative 
to the moving plasma. The speed on the sun-ward side of the shock (i.e., ahead 
of the shock) was v 0n - U = 141 km/sec. This is greater than the. magneto- 
acoustic speed V M1 = 95 km/sec. On the opposite side it was (v — U) = 76 km 
/sec which is less than the magneto- acoustic speed = no km/sec. Thus the 
normal- speed relative the surface shock did decrease, as required for a shock 
and the flow was super magneto- acoustic ahead of the shock and submagneto- acoustic 
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behind it. The component of v tangentto the surface was found to be 414 km/sec 
on side 1 and 416 km/sec on side 2, consistent with the requirement that v lt = 
v 2t . The condition Bj/iij = B 2 /n 2 , from the frozen-field condition, is ’ 
satisfied within the errors for this shock, for Table IV gives 

B^rij = 1.2 ± -2, and B 2 /n 2 = 1.5 + .1 

The change in the' momentum flux was p-v 2 | ^ =8.6 x 10" 10 ergs/cm 2 . This was 
balanced by a change in the pressure 

1 

P s ^ Pi + — — = 8.0 x 10“ 10 ergs/cm 2 

' 1 8 77 

2 

(computed as for the forward perpendicular shock in the preceding section), as 
required by (5.24).) 

Similarly the change in the normal kinetic energy flux, 7.1 x 10 13 (cm/sec) 2 was 
found to agree with the change in the E.HS of (5.25), 5.5 x 10 13 (cm/sec) 2 , within 
the experimental errors. 

Thus, all of the conditions for a reverse perpendicular shock were satisfied. 

5. Fast Shock. Two lands of shocks can propagate at an angle with respect to the 
magnetic field. These correspond to the 2 kinds of magneto acoustic wave, the 
fast wave and the slow wave. Gas enters the fast shock at a speed greater than 
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the speed of the fast wave in the direction of the shock normal and leaves at a 
speed less than this. Thus,, the flow is superfast ahead of the fast shock and 
sub fast behind it, but the flow is super- Alfvenic on both sides of the shock 

J 

(Jeffrey and Taniuti, 1964). The fast shock is the most frequently observed type 
of shock in the solar wind, perhaps because of the similarity of the fast magneto- 
acoustic wave with the acoustic wave in aerodynamics. 

Since the fast shock normal is not parallel to B , B n and B t are both non- zero 
(See Figure 34) Equation shows that B lt is parallel to B 2t at a fast shock; hence 
B x , B 2 andn are eoplanar. Equation 5.10 shows that v - v is parallel to 
B lt - B 2t , and is thus parallel to B. lt and B 2t and eoplanar with n . Then 
v a - v 2 andn are also eoplanar. Summarizing, n , Bj ,B 2 , and v -v are 
eoplanar. This was demonstrated for an anisotropic plasma by Chao (1970) and 
independently by Hudson (1970) and Neubauer (1970). 

The shock normal can be calculated from B, and B„ as follows. Since B, , b 

12 i 2 

and n are eoplanar, B 1 x B ^ is in the plane of the shock. It was shown that 
B - B is also in the shock plane. Thus n is normal to both B, x b„ and 

x 2 12 

Si ' B 2’ 

" = (Ij - B 2 ) X (B, x|)/| (Bj* - B 2 ) X (B, X B 2 ) I (5.26) 

Note that the formula does not give a normal for parallel or perpendicular 
shocks. It can be used to determine n for fast shocks only if the angle co between 

0 
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B i andB 2 , is larger than the uncertainties in Bj andB 2 . Since co is generally less 
than 15° for shocks in the solar wind and since the uncertainties in B 3 and B 2 
are commonly^ 5° or more, the errors in n determined this way may be very 
large, sometimes so large that n is meaningless. Clearly, more accurate normals 
can be obtained by reducing the uncertainties in B x andB 2 . Lepping and 
Argentiero (1971) accomplished this by developing a n least-squares ,T technique 
which gives the best estimate values of B x and 3 2 by requiring that the obser- 
vations of v x f v 2) B p B 2) p x a ndp 2 satisfy (5.1), (5.10), (5.13) with £ ~ 1 and 
L B n l ^ = 0. As an example, consider the shock seen by Ness’ (Pioneer 7) mag- 
netometer on Aug. 29, 1966. Using mean values of & in (5.26) gives an error 
cone for n of 25°. The Lepping- Argentiero method gives an error cone of only 
6°. (The error cone is here defined such that there is a 95% confidence that 
the normal is within it.) Chao (1970). introduced another method for obtaining 
best estimates of j3, an iterative trial and error procedure, which is simpler than 
the Lepping-Argentiero method, but does not always give as much accuracy. The 
accurate determination of shock normals is essential for studies of shock geometry an 
and the interactions of shock with discontinuities. 

Other methods for determining n can be used for shocks that are identified at 
2 or more spacecraft, if one can assume that the shock surface is plane over 
distances equal to the spacecraft separations. Sonett et al. (1964) used such a 
method to determine n for the shock on Oct. 7, 1962. Chao (-1970) gives the 
formula n = [b 3 x (v 2 -v 3 )] x U/| [B x x (v 2 - v 1 )]-x U | where U =R /T - 
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(P 2 v 2 - p 2 v 2 )/{p 2 ~ p 1 ) s Rand T being the vector displacement and time delay 
between the spacecraft. Ogilvie and Burlaga (1969, erratum, 1970) introduced a 
method for calculating n using simultaneous measurements from 3 spacecraft, 
the four parameters n,U were determined from the 3 arrival times and the 

condition that B - B must lie in the shock plane. 

2 1 

Greenstadtet al. (1970) used onset times from 3 spacecraft and certain limits on 
the shock speed to obtain limits on n for the 7 July 1966 shock. These limits 
were narrowed using the time of the associated sudden commencement. Their 
basic result, that the shock normal pointed 65° to 70° below the ecliptic, is 
surprising. It should be noted that the 3 spacecraft which they used were all 
close to the ecliptic plane, so the most important component of n , n , was the 
least accurately determined. If one uses the intercept of shock surface with the 
ecliptic (T = (.276, .961, 0)) which was determined by Greenstadt et al. (which 
is quite accurately determined by their method since the spacecraft were widely 
separated in the ecliptic plane) and if one uses the value B i - B ? = (-2.0, 6.4, 4.4) 
which they give, then the method of Ogilive and Burlaga (1969) gives n = T x 
(S x -B 2 )/|T X (B - B 2 ) [ - (.96 - .275, - .04), which implies an inclination 
$ n = - 40°. This smaller value of Q (i.e., shock normal closex- to the ecliptic 
plane) is more consistent with other values that have been published (Chao, 1970, 
Ogilvie and Burlaga, 1969), and agrees closely with the value 6 n = - 37° obtained 
from (5,26) by Greenstadt et al. 
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Another method of determining shock normals was given by Burlaga (1970) who 
used onset times from 3 spacecraft, the measured solar wind speed, and (5.1) to 
calculate n and U by means of an iteration scheme. In his equations (13) and (14), 
tan 6 should be replaced by tan 9 sec <f > . 

Given the shock orientation and a complete set of measurements, one can easily 
determine whether the momentum flux and energy flux equations are satisfied. 
However, in practice one cannot measure all of the quantities in these equa- 
tions. The magnetic field is known more precisely than the plasma parameters, 
so it is convenient to replace these equations by expressions for the jumps n 2 /n 1 , 
T,| 2 /T 1U , T i2 /T 1 as functions of the magnetic field parameters and. the upstream 
Alfven Mach number. Using frozen-field conditions (5.3), the transverse 
momentum flux equation (5.5) can be written 


fix 

Pi 



B 




B 


it 


B. 


2 1 . 


3 


2t 


(5.27) 


This is equivalent to (5.6). Chao (1970) has used a form of this equation to study 
the variation of nj/^ versus for various values of the other parameters. 
When ^ = ^2=1 (5.27) reduces to an equation derived by Wilkerson (1969). 

For a gas consisting of protons and electrons the normal momentum flux 
equation (5 .4) can be written 



(5.28) 
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where 

, />i(V ln --U ) 2 

f ° "i k (1* + IT,) 

B? 

and «-i - i 

8 7r n 2 k (Tj^ + Tjp 

g jl 

~~ 4 tt n x k (T+ + Ti -) 

This also reduces to an equation derived by Wilkerson (1969) when ^ 2 ~ 1. 

Now let us turn briefly to the observations of fast shocks . Numerous fast for- 
ward shocks have been identified in the solar wind, the first by Sonett et al. 
(1964). A list of fast shocks and the corresponding references may be found in 
Hundhausen (1970). We shall give just one example, the shock of Aug. 29, 1966. 

Plasma data from the MIT probe on Pioneer 7 and simultaneous magnetic field 
from Ness’ magnetometer are summarized in Figure 35 from Chao (1970) to- 
gether with Chao’s "best fit" to the Rankine-Hugoniot equations with = £ 2 . 
Table IV gives the somewhat better "best fit” values of Lepping and Argenterio 
(1970). Clearly, the R-H equations are approximately satisfied. Lepping and 
Argentiero, using the method described above, found that the shock normal was 
4 > = 197°, 0 - 19°, within an error cone of 6°. Using this value for h and the 

n n. 

time 48.5 min that it took for the shock to move from Explorer 33 to Pioneer 7, 
which were separated by the vector R = (1.29, .59, .16) x 10 6 km/see, (See 
Figure 36). The shock speed was found to be 471 km/sec. This compares well 
with the value 467 km/sec obtained by Lepping and Argentiero from the Rankine 


Hugoniot equations. 
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The observation of a reverse fast shock was recently reported by Chao and 


Binsack (1971), The earth's bow shock is a stationary reverse fast shock which 
propagates toward the sun at a speed just equal to the solar wind speed. It has 
been extensively studied, and the. results are summarized in recent reviews by 
Spreiter and Alksne (1969) and by Wolfe and Intriligator (1970). 

n a = n 2 shocks. The above discussion assumes that n t ^ n 2 for the fast shocks. 
Formally, one can set n = n and still have v ^ 0. in this case, v = v . 

Now if V A1 / v; 2 , then (5.7) shows that B lt and B 2t are colinear, as for a fast 
shock. Solving (5.7) for V n - U gives 


(v „U) 


In 


2 _ V a? 8 i, - Vii B 2 , 

B I. - B 2. 


Introducing the definition of V ' and <f , gives 


(5.29) 


n fC B it - B 2t ) 5 i B 2 + 477 B 2t CP|! 1 -P l2 ) B2 -4rrB 

Clearly, B o for these discontinuities. 

Shocks of this sort were predicted by Ivanov (1970) who called' them "rotational 
discontinuities" because in the limit of zero anisotropy the shock speed goes to 
V A . This name is somewhat misleading, however, since an essential property 
is the colinearity of B lt and B 2t , which makes it more like a fast shock than 



it 2 


(p 


111 
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a rotational discontinuity. (Ordinary rotational discontinuities derived their name 
from the characteristic that B t rotates across the plane of the discontinuity so 
that B 2t and B lt .are generally not parallel.) It is true that 



for both Ivanov-’s discontinuities and Alfven shocks, but this follows directly 
from (5.5) and the assumption that p l = p 2 . Since B and B 2fc are colinear, the 
coplanarity theorem applies to n x = n 2 discontinuities (but not to Alfven 
shocks), and the normal is given by (5.26). 

Ivanov (1970) suggested that 10 out of 11 of the large velocity discontinuities 

in Burlaga (1969b) are shocks with n =n . However, additional data presented 

$ 

by Burlaga (1971) show that this cannot be correct, since in the cases where 
the orientation of the surfaces can be determined directly, it is not that which 
is predicted by Ivanov. Further evidence in Burlaga and Chao (1971) shows that 
discontinuous changes in B are accompanied by density changes at least 95% of 
the time, which indicates that n 1 = n 2 shocks are seldom, if ever, present. 

6. Slow Shocks. The slow shock propagates at an angle with respect 
to B (see Figure 37) and corresponds to the slow magneto-acoustic wave mode. 
Gas enters at a super slow speed and leaves at a subslow speed. The flow is 
sub-alfvenic on both sides of the shock. An important and distinctive character- 
istic of a slow shock is that the magnetic field intensity is lower behind the shock 
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than ahead of it; magnetic energy is lost, presumably by. conversion to thermal 
energy or possibly to wave energy. The fractional change in 33^ , ~ [B^j/B^ , is 
a measure of the shock strength. 

The Rankine-Hugoniot equations for slow shocks are essentially the same as for 
fast shocks. 

Slow shocks were first identified in the solar wind by Chao and Olbert (1970) us- 
ing Mariner 5 plasma and magnetic field data. The observations of the 2 slow 
shocks discussed by Chao and Olbert are shown in Figure 38. It is clear that 
these slow shocks are not as thin as fast shocks and the fluctuations near these 
slow shocks are larger than is usually observed at fast shocks. Nevertheless, 
the signature of each of these 2 discontinuities is clearly that of a slow shock 
(an increase in n, v and T and a decrease in B). The dashed lines in Figure 38 
are "best- fits" to the Rankine-Hugoniot equations, obtained by the trial and 
error procedure of Chao (1970). These are not unique fits, but they do show 
that there is a solution of the Rankine -Hugoniot equations for an isotropic 
plasma which is consistent with the observations. The condition that the flow 
must be subalfvemc on both sides of a slow shock was satisfied for the 2 
shocks in Chao and Olbert (1970). The condition that flow be superslow ahead 
and sub slow behind the shocks was not verified directly because the electron 
temperature was not known; however, using the value of determined by 

their fit to the Rankine-Hugoniot equations, Chao and Olbert found M j =1.5 
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before and .7 after the shock on July2Q, 1967, and M slow =1.8 before and .6 
after the shock on August 30, 1967. 

Another slow shock was identified in Pioneer 6 data by Burlaga and Chao (1971). 
This is shown in Figure 39 together with fits to the R-H equations obtained by 
Chao's method. The speed normal (on RTN coordinates) and Mach numbers for 
this shock are shown in Table V together with the same parameters for the 
shocks in Chao and Olbert (1970). 

A reverse slow shock was also identified in the Pioneer 6 data by Burlaga and 
Chao (1970). The data are shown in Figure 40 together with a fit to the Rankine 
Hugoniot equations using the method of Chao. The speed, orientation and Mach 
numbers for this shock, which occurred at 0901 UT on Jan. 19, 1966, are shown 
in Table V. There was only a small drop in the flow speed across the shock, 
from 29 km/sec to 23 km/sec, but there was a significant drop in the magnetic 
field intensity, from 5.2y to 4.2y. The loss of magnetic energy is comparable to 
the loss in flow energy. 

The existence of slow shocks in the solar wind is strong indirect evidence for 
the existence of the slow magneto- acoustic wave and the ion-acoustic wave in the 
solar wind. 


6. Alfven Shock. Strictly speaking there is no Alfven shock correspond- 
ing to the Alfven wave in the same sense that fast and slow shock correspond to 
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Table 5 


Slow Shocks 


Date 

July 20, 1967 

Forward 
Aug. 30, 1967 

Jan. 20, 1966 

Reverse 
Jan. 19, 1966 

h 

(.38, -.80, .47) 

(.01, -.09, -.99) 

(.75, .30, .60) 

(.85, .01, .54) 

V (km/sec) 

130 

2 

323 

251 

M a,i 

0.9 

0.9 

0.9 

0.9 

M a,2 

0.7 

0.6 

oo 

© 

0.8 

^■slow, 1 

1.5 

1.8 

1.3 

1.2 

2 

0.7 

0.6 

0.8 

0.8 


fast and slow magneto acoustic wave. This point is clear from our earlier dis- 
cussion of parallel shocks. It can probably be traced to the fact that Alfven 

/ 

waves are not compressive. There is, however, another shock (G ^ 0) which is 
unique to magnetohydrodynamics; it is called an Alfven shock. This is illus- 
trated in Figure 41. The flow is alfvenic on both sides of this shock (Mj^ = 

M ' — 1), so the discontinuity propagates at the Alfven speed, essentially like a 

kink in the magnetic field. For an Alfven shock B n 0, and (5.7) shows that 
B lt and B need not be parallel (see Figure 41). In fact, the Alfven shock is 
often called a rotational discontinuity because B 2t is rotated in the shock sur- 
face with respect to B . Clearly, b 1 , B 2 and n are not coplanar, in distinction 
to the fast and slowshoeks. In this case, B 2 x B 2 = B n x (B 2t - B lt ) + B lt x 
B„. is not in the shock surface, so the Alfven shock normal is not gi'/er by 26 ). 

/V 1 1 
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Equation (5.13) shows that the velocity change and magnetic field change are re- 
lated by 


-Si ^1 _ S2 ^2 ^ 

,\/4 if p x \/ 4 7 t p 2 J 


(5.31) 


Since 


B* 


02 = p 2 (v _ u)2 = p 2 v'2 = p £ 


4 77 


the normal momentum flux equation (5.12) gives 




= 0 


which is the same as the condition derived for a tangential discontinuity, 
energy flux equation (5.14) becomes 


(5.32) 


The 


e + Z](P|| + PjV2 + B 2 /8 
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= 0 


(5.33) 


L_ -n 

If the anisotropy does not change across an Alfven shock, £ - £ 2 ’ ^ en the con- 
servation of mass flux (5.1) implies that p 1 = p 2 . Then (5.32) and (5.33) give 
£ i = e 2 which implies that p 2 = p.. Putting this into (5.32) gives B. = 8„ . 

Thus, if = f 2 , then = B 2 and p x - p 2 - However, if ^ ^ £ 2 , B 1 ^ B 2 ; this 
important fact was first noted by Hudson (1970). 


Evidence for Alfven shocks was found by Belcher and Davis (1971) in the 

\ 

Mariner 5 plasma data of Bridge, Lazarus, and Snyder and the magnetic 
field data of Coleman, Jones and Smith. This is shown in Figure 42, 
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/here it is seen that (5.31) with g x = £ 2 = 1 is satisfied for at least 2 discon- 
inuities. There was no appreciable change in n, as required for an Alfven 
hock with f x ~ No temperature measurements were given, but presumably 
ET X = 2 T 2 . It was found that B = B 2 , but this result is not independent of an 
ts sumption which was made to obtain the spacecraft field. 

Alfven shocks may be pictured as propagating "kinks" in the interplanetary 
nagnetic field. Parker (1963) was the first to suggest that they are present 
ind indicated that they might be important as cosmic ray scattering centers. 
Juenbyet al. (1970) have given this idea some experimental support. No 
neasurement of the speed of an Alfven shock has been made. 

Classification of Shocks, Here we simply note how the types of shocks 
liscussed above can be viewed as different solutions of (5.7), which can be 
vritten 


= P 2 Kl( 1-M$ 

Parallel shock: B = B = o 

it 2t u 

Perpendicular shock: v ' 2 - v ' 2 - 0 
Alfven shock: M' = 1; M' = 1 
Fast Shock: > 1; M^ 2 > 1 
(n x = n 2 shock. p l =p 2 ' ) 

Slow shock: < 1; < 1 


(5.34) 
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Other solutions are formally possible: = 1, £ 1 (switch-off shock). 

£ 1, M a ' 2 = 1 (switch-on shock); > 1, M A ' 2 < 1; M; 2 < 1, > 1. 

i 

But these are not "evolutionary", and they are not expected to be observed in 
the solar wind. 

YI. Summary and Discussions of Problems for Future Study 
We have considered the applicability of hydromagnetic theory to the study of the 
microstructure of the solar wind. In particular, we have presented the theory of 
hydromagnetic waves and discontinuities which is appropriate for the solar- 
wind, and we have reviewed the experimental evidence for the various waves, 
discontinuities and some of the instabilities which are predicted by this theory. 

Nearly all of the discontinuities given by the theory have been shown to exist in 
the solar wind. These include tangential discontinuities, forward and reverse 
fast and slow shocks, perpendicular shocks, and Alfven shocks. Parallel shocks 
have not been found; their existence would be of special interest, since it would 
be clear evidence for the existence of a longitudinal wave mode and would show 
in the simplest form the way that the anisotropy changes across a shock. Con- 
tact discontinuities have also not been found; their discovery would tell us some- 
thing about the diffusion of particles along B. The origin of all types of discon- 
tinuities (except some of the fast shocks, which are caused by solar- flares) is 
one of the outstanding problems . The mutual interaction of discontinuities and 
the intersections of discontinuities is justbepnning to be studied. The structure 
of discontinuities is an important subject which can only be studied ; by multiple 
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spacecraft observations with high time resolution instruments; it is not clear 
to what extent this problem can be treated by the MHD approximation. 

The study of waves is in a somewhat less satisfactory state. Alfven waves have 
been identified with some certainty, but their propagation speed and the dispersion 

relation have not been measured. It is not clear whether or not an acoustic wave 
can propagate along B . The existence of slow shocks in the solar wind and the 
theory of ion acoustic waves in a collisionless plasma indicate that acoustic 
waves do exist in the solar wind, but they have not been measured directly and 
there are theories which suggest that particle motions parallel to B are not 
strongly coupled, in which case acoustic waves could not propagate. Sinusoidal, 
compressional oscillations have been observed, which may be magneto- acoustic 
waves, but more precise density and velocity measurements will be needed to 
positively identify such waves. Evidence for the fast magneto-acoustic mode, 
propagating anisotropically as predicted by MHD and particle orbit theory, has 
been found in the lunar Mach cone. However, it has not been established that 
the corresponding slow mode can propagate. Future studies of flow past some 
of the outer planets will be especially significant in this regard if the solar wind 
becomes transonic or subsonic there and retains its characteristics as a hydro- 
magnetic fluid. The firehose instability has often been discussed in reference 
to the solar wind, but its existence has not been firmly established. It is clear, 
however, that the instability is not as universal as some authors have suggested. 
More work needs to be done in this area. 
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The experimental results which we have reviewed show without doubt that 
hydromagnetic theory is applicable to the solar wind. Many of the consequences 
of this theory remain to be worked out. Nonlinear theories of waves, instabilities, 

and -interaction of discontinuities and waves are particularly needed. One would 
also hope for future developments and solar wind applications of the MHD theory 

of turbulence, boundary layers, and flows past obstacles. Experimentally, the 
need is for higher time resolution, greater precision and simultaneous measure- 
ments of B and the distribution functions of protons, a’s and electrons from at 
*\> 

least 2 space probes. 
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Figure 1. This defines the scales in terms of characteristic times, 
frequencies, and lengths. Some characteristic features seen on the 
different scales are shown. E is the energy of a proton whose gyro 
radius is the scale length in a 5 y field. 
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MAGNETIC FIELD, 



Figure 3. Alfven wave. The perturbations 
Sv and SB are perpendicular to k and B 

~ <-o -v 



TIME, 36.96 SEC UNITS 

Figure 4. Observations of an Alfven wave. This shows a nearly 
periodic wave form with a strong correlation between the radial 
component of the bulk speed (Vp) and the radial component of 

B(Bx). 
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Figure 5. Alfven wave. The observed velocity and that computed from B assuming an Alfven wave for the wave 
in Fig. 4 are in good agreement. The corrected density is constant, as required for an Alfven wave. 
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Figure 6. Observations of non-linear, a periodic 
Alfven- waves by Belcher et dl. These waves 
have the expected proportionality constant be- 
tween Sv^and SB as indicated by the top distri- 
bution. The bottom di stribution'shows that they 
are transverse. 
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Figure 7. Magnetoacoustic Wave. This propagates normal 
to B. Sv is normal to B and '§ B is parallel to B. 
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bly a magnetoacoustic wave. (See text). 














Figure 13. Evidence for coupling between fast and Alfven waves. This shows that the anisotropy of the 
magnetic fluctuations/measured by 7 ], is smallest when the amplitude of the transverse and longitudinal 
fluctuations, measured by A B, to is largest. 
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Figure 14- Test for firehose instability. The instability condition [ 3 ^ — (3 > 2 
was not satisfied for the period shown here. 
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Figure 18. Construction of Mach- cones. There 
is one cone for the fast mode and one for the 
slow mode. Theheavy line is the fluid velocity. 
The '“cones” are constructed by drawing tan- 
gents to the wave normal surfaces from this vec- 
tor. 
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Figure 20. Observations of three tangential discontinuities. The discontinuities are indicated by the dashed lines. 
The proportionality between and n throughout the interval indicates that the pressure is constant and does not 
change across the discontinuities. 







Figure 21. A tangential discontinuity observed at 2 widely 
by 57. 5rhin, thedelay predicted for the orientation computed 






Figure 22. Multi satellite observations of discontinuities. Each panel represents one 
discontinuity. Each discontinuity (represented by aline segment) was seen at 3 space- 
craft. The observed time delays are consistent with the orientations computed assum- 
ing that the discontinuities are tangential. 
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Figure 23 (b). For a tangential discontinuity the orientation may be 
computed from B and independently from v Observed discontinuities 
are shown by dots. 
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RELATIVE SCALE 



r( sec): 

Figure 24. ‘'Thickness" distribution (f D (t)) of discontinuities'. TheactuabtKickness 
is on the order of the solar wind speed times X Most discontinuities* tend to be very 
thin, and are .distinct from waves. 
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Figure 25. Current sheets associated with discontinuities. The magnetic field vectors 
tends to rotate in a plane as shown atthe top of each panel. There are fluctuations out 
of this plane, as indicated at the bottom of each panel, but these are not statistically 
significant. 
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Figure 26. Laminar and “turbulent” 'current sheets. The top panel 
shows2 current sheets on which the field rotated smoothly from its 
position on one side of the discontinuity to that on the other side. 
The bottom panel shows an irregular transition, suggesting a “tur- 
bulent” or unstable current sheet. 
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Figure 28. Annihilation at D sheets. The annihilation 
hypothesis implies that B min /B = cos {co/2) where co is 
the change in the field direction. This seems to be the 
case for the D sheets in the top of this figure. 















‘Figure-31.' Perpendicular ;sh'ock. 
B i s .perpendicuiar*to- r rr. 



Figure 32. Observation of a perpendicular shock. The shock orientation and speed 
shown here were computed from the arrivaltimes atthe 3 spacecraft which are shown. 
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Figure 34'. Fast Shock. A special 
case-as shown for simplicity. 
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Figure 35. Observations of 2 fast shocks: The -solid horizontal segments are 
are averages, the dashed segments are Chao’s “best-fit” values. 
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Figure 36. Observations of a fast shock by 2 spacecraft. The normals sHoton here 
were computed by Lepping and Ar'gentiero. 
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Figure 37. Slow shock. A 
special case is shown. In 
general one can add a com- 
ponent of V in the shock 
plane. 
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Figure 38. Observations of 2 slow shocks. Observed (solid lines) and 
calculated values (dashed tines) are in good agreement. 
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Figure 41. Alfven shock. 
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